Generation of electron spin polarization in disordered organic semiconductors 
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The generation mechanisms of electron spin polarization (ESP) of charge carriers (electrons and 
holes, called "doublets") in doublet-doublet recombination and triplet-doublet quenching in disor- 
dered organic semiconductors are analyzed in detail. The ESP is assumed to result from quantum 
transitions between the states of the spin Hamiltonian of the pair of interacting particles. The 
value of the ESP is essentially determined by the mechanism of relative motion of particles. In 
our work we have considered the cage and free diffusion models. The effect of possible attractive 
spin-independent interactions between particles is also analyzed. Estimation with obtained formulas 
shows that the proposed mechanisms can lead to a fairly strong ESP much larger than the thermal 
one (at room temperatures). 

PACS numbers: 73.50.-h, 73.43.Qt, 75.47.-m, 72.25.Dc 



I. INTRODUCTION 

Electron spin polarization (ESP) of charge carriers 
(electrons and holes) in solid-state systems is the im- 
portant phenomenon, which is observed in ferromagnetic 
metals and inorganic semiconductors and result from 
strong magnetic interactions in these solids*^ The ESP 
effect is very interesting from general physical point of 
view and important for applications**— 

As for organic semiconductors (OSCs), magnetic inter- 
actions in them are typically fairly weak and no strong 
ESP is expected. The ESP can, nevertheless, be in- 
jected from ferromagnetic solids^ In many experiments 
the OSCs are applied as spacers (between ferromagnetic 
leads) , in which one can manipulate the ESP by applying 
electric and magnetic fields The dynamics of electron 
spins in OSCs is studied quite thoroughly^ but the ma- 
jority of works concern the analysis of spin evolution of 
polarized charge carriers injected into OSCs from polar- 
izing leads (see, for example, refs. [4-10], as well as [2,3] 
and references therein). 

This work concerns the discussion of the mechanisms of 
ESP generation directly in OSCs (in the absence of spin 
injection). The possibility of this effect is usually ignored 
because of above-mentioned weakness of magnetic inter- 
actions in these materials. The analysis shows, however, 
that some spin selective processes in OSCs can result in 
fairly strong ESP (much larger than the thermal one). 
Such processes are similar to those, which determine the 
magnetic field and spin effects (MFEs) in spin selective 
reactions of paramagnetic particles [usually doublets (D) 
and triplets (T)]^M 

In our investigation we will consider two types of pro- 
cesses, doublet-doublet (D — D) recombinatio n 14 ! 15 and 
triplet-doublet (T-D) quenching in which the ESP 
generation is expected to be very efficient. The ESP 
results from quantum transitions between the terms of 
the spin Hamiltonian of the pair of particles (D — D or 
T—D), assisted by their stochastic relative motion. The 
terms are determined by the Zeeman and spin exchange 



interactions and transitions are induced by the hyper- 
fine (HFI), anisotropic part of the Zeeman (AZI), or zero 
field splitting (ZFS) interactions, depending on the pro- 
cess considered. The strong dependence of the exchange 
interaction on the interparticle distance r leads to the 
significant localization of transitions in space and essen- 
tial dependence of the ESP generation efficiency on the 
characteristic properties of relative motion. 

Noteworthy is that the theory of diffusion assisted ESP 
generation in liquid phase reactions has already been 
developed*^ - — It is based on the stochastic Liouville 
equation (SLE) approach ) 22 ' 23 which allows one to rig- 
orously describe the generation of MFEs, (including the 
ESP). As applied to similar processes in OSCs, however, 
this SLE-based theory needs certain modifications, due 
to some specific features of interaction and relative mo- 
tion of paramagnetic particles in these materials. 

The effect of these specific features on the ESP in OSCs 
is considered within two models of relative motion of par- 
ticles: cage and free diffusion. Both models are shown 
to predict similar expressions for the ESP both in gem- 
inate and bulk processes. The proposed theory enables 
one to take into account possible spin- independent (i.e. 
Coulomb) interaction U(r) between reacting particles, 
which can strongly affect the MFEsj 24 ' 25 Simple estima- 
tions demonstrate that the attractive interaction U(r) 
results in significant increase of the ESP. The evaluation 
of the ESP for realistic values of parameters of the the- 
ory shows that the considered mechanisms can be very 
efficient leading to a quite large ESP in OSCs. 



II. GENERAL FORMULATION 

The mechanisms of the ESP generation in diffusion as- 
sisted liquid pase spin- selective reactions are thoroughly 
investigated in the theory of MFEsfiI~— In our work we 
will briefly analyze most important results, but will be 
mainly concerned with extending and adapting them to 
the processes under study. 

We consider two types of spin-selective processes: 
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the doublet-doublet (D — D) recombination and triplet- 
doublet (T—D) quenching, in which doublet particles are 
associated with electrons or holes (D = e,h). The most 
part of the consideration concerns geminate reactions, 
though bulk processes are also discussed. 

The kinetic scheme the geminate variant of both pro- 
cesses can be represented as 

X a + X h ^- [X a ...X b ] ^X R , (X = D,T). (2.1) 

Here X a and Xb denote D- or T-particles involved in the 
reaction. The reaction kinetics is essentially controlled 
by the rate w d of dissociation of the coupled state (cage) 
[X a ■ ■ ■ X ] and spin selective reaction rate K R . 

The kinetics of ESP generation in these processes is 
described with the spin density matrix pit) of pairs of 
paramagnetic particles X a and Xb- Naturally, pit) de- 
pends on the intcrparticle coordinate r. For simplicity, 
we assume that the precess (|2.1[) is isotropic, i.e. the in- 
terparticle interactions and initial condition depend on 
the distance r = |r| between particles. In this case p 
depends only on r: p(r, t) = p(r, t) 

The spin evolution of the pairs is governed by the spin 
Hamiltonian H(r,t), which can, in general, be time de- 
pendent. The explicit form of H(r,t) will be specified 
below for both considered processes. 

In the above- formulated assumptions the space/time 
evolution of the pair density matrix p(r, t) is determined 
by the SLE (h = 1) 



p=-i[H,p]-(K R +L) P) 



(2.2) 



with [H, p) = Hp — pH. In this SLE the operator L 
describes the relative motion of particles, which is treated 
as hopping migration and is characterized by the average 
jump length A and the rate of jumps Wh- The term K R p 
represents the effect spin selective reaction: 



Kp = W{p r p+pp r ) 



(2.3) 



where P is the operator of projection onto the reactive 
states of the pair and k r is the reaction rate. In general, 
k r = K R (r) depends on the distance r and this depen- 
dence n R (r) is strongly localized at r small distances. 

The ESP essentially depends on the initial spin and 
spatial state. In general, the initial density matrix can 
be written as 



pit = 0) = p (r)p Q , 



(2.4) 



where p (r) is the initial spatial distribution (see below) 
and p is the initial spin density matrix. 

The observable under study is the ESP V(t — » oo) of 
-D-particles. For brevity of notations, in what follow we 
will mainly analyzed the absolute value of the ESP P. 
As for the ESP sign, it can easily be obtained, as it is 
demonstrated in Sec. V. 

In our work we consider the most important mech- 
anisms of the ESP generation in two above-mentioned 



processes and within cage and free diffusion models of 
relative motion of particles: 

1) Cage model. In the cage model the considered pro- 
cess is represented by the kinetic scheme (|2.1[) , in which 
w d is the rate of monomolecular dissociation of the cage 
and K R is the rate operator of spin-selective reaction, 
independent of the distance r (n R is independent of r). 

The operator L, essentially controlling the kinetics of 
the process [see eq. (|2.3|) ]. is defined by 



Lp = w d p. 



(2.5) 



Detailed analysis shows 18 that non-adiabatic transi- 
tions between terms of the Hamiltonian H, which deter- 
mine the ESP generation, are localized near the distance 
d of closest approach (of reacting particles) in the narrow 



region of width Sr 



■C d, where a is the size of 



the exchange interaction (see below). The distance d is, 
in turn, expected to be of order of the hopping length A. 
In such a case, it is quite reasonable to assume, that the 
cage size r c ~ A (~ d), and associate the cage dissocia- 
tion rate w d with the rate Wh of jumps: w d ~ w h . 

For pairs in these cages (of volume V c ~ A 3 ) the rel- 
ative X a — Xb coordinate r is suggested to be randomly 
(homogeneously) distributed. Note that for the homo- 
geneous distribution over r the initial condition for the 
SLE can be written in the form (|2.4I) : pit = 0) = p a po 
with p = V~ 1 . 

In the cage model, for example, the ESP of Z?-particles, 
escaped from the cage, can be evaluated by formula 



P = 2w 



dt(TT[S DzP (r,t)}) reVc 
2w d \(Tr[S at p(r,€ = 0)]) v \, 



(2.' 



where S D is the z-projection (along the magnetic field 



B) of the £>-particle spin S D , p(e) 
Laplace transform of pit) and 



J Q dt e et p(t) is the 



(P)r 



v: 



d 3 r p(r) 



(2.7) 



is the average over the homogeneous distribution in the 
cage (of volume V c ). 

2) Free diffusion model. In the free diffusion model the 
considered processes can, formally, be expressed by the 
same scheme (|2.1[) , in which, however, the reaction and 
dissociation kinetics is non-exponential and controlled by 
relative diffusion of reacting particles. In this model 

Lp = -D p r~ x V^{rp) and k r = K r S(r - d), (2.8) 



where D p = D x 



D Xb , (D p ~ X 2 Wh), is the coefficient 



of relative diffusion of particles X a and Xb, V r = d/dr, 
and k R is defined by k R = ^K r (P p + pP ) [eq. (|2.3|) ]. In 
the free diffusion model the SLE should be solved with 
the reflective inner boundary condition at the distance 
of closest approach di^ V r p(r,t)\ r= d = 0, and the outer 
condition p(r oo,i) — > 0. 
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The ESP of Z?-particle is expressed by formul a 18 ' 19 



V = 2 



d 3 rTv[S D p(r,t^ oo)] 



(2.9) 



The kinetics of ESP generation in the free diffusion 
model has already been analyzed in detail for both D—D 
and T — D processes^— that is why in what follows we 
will restrict ourselves to brief summary of some results 
of this analysis interesting for our work. 




III. D — D RECOMBINATION 



The D — D-recombination is described by the kinetic 
scheme (|2.8|) with X a = D a and Xb — Db, in which K R 
is the operator of reaction in the singlet (S) state of the 
D a -Db pair, given by eq. with P R = P g = \S){S\. 

The most important mechanism of ESP generation in 
this type of reactions is known as 5T_-mechanismi 14 i 15 
Here we will consider the specific features of this mecha- 
nism in organic semiconductors. 

A. Interactions 

The spin Hamiltonian of the D a — Db pair is assumed 
to be represented as a sum 

H(r,t) = H (r) + V(t), (3.1) 

where Hq(t) and V(t) are the steady state and randomly 
distributed (and may be fluctuating) parts, respectively. 
The steady state part of the Hamiltonian 

H (r) = H z + H ex (r) (3.2) 

is diagonal in the basis of the total electron spin S = 
S a + Sf, (called ST-basis) and consists of the Zeeman 
interaction (in the external magnetic field B) 

H z = n (Sa z +S b J, (3.3) 

in which f2o — gPB is the Zeeman splitting (we assume 
that the isotropic parts of ^-factors of both particles are 
the same: g a = g~b = g = 2), and the exchange interaction 

H ex = J(r)(i + 2S a S b ) (3.4) 

exponentially decreasing with r: 

J{r) = J cxp[~a{r - d)}, (J > 0), (3.5) 

where d is the distance of closest approach: d 3> or 1 . 

In the ST_ -mechanism the ESP is assumed to re- 
sult from quantum transitions between U s (r) and U T (r) 

terms of H (r) near the 531 -crossing region at r = r t 
(see Fig. 1). These transitions are caused by the fluctu- 
ating part V(t) of the Hamiltonian (non-diagonal in the 
ST-basis), which is a sum of two terms: 

V(t)=V a (t) + V z (t), (3.6) 



FIG. 1. The schematic picture of terms U M (r), (M = 
S,T a ,T ± ), of the spin Hamiltonian H (r) [eq. P^]) ] of the 
D — D pair (for J > 0), in which r t is the coordinate of the 
region of efficient quantum transitions. 

The first term V H (t) represents the hyperfine interac- 
tion (HFI), which in the considered realistic case of a 
large number of nuclei can quite accurately be treated as 
the interaction of electron spins S a and S;, with (classi- 
cal) random magnetic fields B Q and B b , respectively: 

V H {t) = g)8[S B (*) + S b B b (t)]. (3.7) 

Distributions of fields are B a and Bf, isotropic, and 
Gaussian with (B„) = and mean squares (B^) ~ 
J2j luj (1 + I^ j ) a t j i iy — °j b), determined by hyperfine 
coupling constants a Vj . 26 Hopping of D a and Db par- 
ticles results in sudden changing of nuclear magnetic 
fields, which can be considered as stochastic vectors 
B„(i), (y — a,b). For simplicity, the correlation func- 
tions of projections B Vq (t) 1 (q — x,y,z) are taken in 
the form {B Vq (t)B Vql (0)} - 6 qq >(Bl) exp(-w H J) (so that 

(V H (t)V H (0))~e- J «S). 

The second term V z (t) in eq. (|3.6|) is the anisotropic 
part of the Zeeman interaction (AZI): 

V z (t)=p{S a g , a (t)B + S b g' b (t)B}, (3.8) 

where 

g' v = 9u - 9v with g v = Tr(g v ), (v = a, 6). (3.9) 

The representation (|3.8p implies that Tr(gl) = 0. As 
in the case of the HFI the time dependence g' v (t) re- 
sults from Z?-particle hopping. We will use the simple 
form of the correlation functional (g' u f (t)g' v / (0)) ~ 

e - w z»\ (q q'), (which means that (V z (t)V z (0)) ~ 



B. ESP generation kinetics 

In order to evaluate the ESP V we need to specify the 
initial spin state p [see eq. (|2.4p ]. In our analysis, for 
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certainty, we consider the equilibrium (isotropic) state: 
Pa = j(Ps+P T ), (3-10) 

where P g = \S)(S\ and P T = £ m=0 ± \ T J( T J- Notc ' 
that for other initial spin states the ESP is, actually, 
proportional to that for p (|3.10[) i2£ 

1. Cage model 

We have already noted that quantum transitions, re- 
sulting in the ESP, are localized in the narrow region of 
width Sr ~ aT 1 <C d (recall that d ~ A). In such a case 
the ESP generation can properly be described by treat- 
ing spin evolution of D—D pairs in all positions of them 
within this region as cages, whose dissociation rate is as- 
sociated with the hopping rate w h (see Sec. II): w d <w h . 

The interaction V in thus introduced cages should be 
considered as time independent though randomly dis- 
tributed, since the fluctuations of V(t) are formed just as 
a result of hopping motion, i.e. at times t>T h = w7 . 

Under the assumption of homogeneous distribution of 
D—D coordinate r in the cage (see Sec. II) and with the 
use of eq. (I2.6P one obtains the following formula for the 
ESP V, generated in D — D recombination, 



\V_ - V, 



l<P- + (r)>, 6 



(3.11) 



where V ± = w d {p ± (r)) rev are the contributions of T±<r^ 
S transitions, in which p±(r) = (T ± |p(r, e = 0)|T ± ), and 
p_ + (r) = w d [p_{r) -p + {r)]. 

Solution of the SLE in the second order in the ran- 
domly distributed interaction V yields formula 



p_ + (r)=p(\V ± \ 2 ) 



1 



1 



+ fi2(r) w 2 c +n 2 {r) 



\(S\V\T + )\ = \(S\V\T_)\, fi ± (r) 



(3.12) 



fio ± 



with | V ± 

2J(r), w c — w d + ^k r being the S— T-dephasing rate in 
the cage, and p = \(k r w c ) /[w d (w d + k r )] 

Substituting this formula into eq. p. lip one can obtain 
the following expression for the SPE V, corresponding to 
the initial equilibrium spin state (|3.10j) : 



P=p \<p u -^> c /(l + w c 2 )], 



(3.13) 



It also is worth noting that the ESP for any other initial 
spin state is closely related to that (13.131) for the equilib- 
rium state p E (|3.10j) . For example, for the triplet initial 
state Po = p T = \P T the ESP V(p T ) = (A/3)V(p E )/p Rc . 

Formula (|3.13p is fairly cumbersome. Fortunately, in 
applications it is quite sufficient to use simpler limiting 
expressions at large Jq > fio, w c and small Jq < Qq,w c : 



V = V S = (irp c ) 



V = V,„ 



8 Pc 



for J > fio, w c , (3.15) 



for J < fi , u> c . (3.16) 



It is easily seen that the considered two cases Jo > fio 
and Jo < fio correspond r t > d and r t < d, respectively, 
where r t is the distance of ST_ crossing (see Fig. 1). 



2. Free diffusion model 

Simple formulas for the ESP have also been derived in 
the free diffusion modelJ^r— Naturally, the ESP, gener- 
ated in diffusion controlled D—D recombination, depends 
(but weakly) on the initial distance r i . Below, for cer- 
tainty, we present the formulas in the realistic case r i « d. 

The most interesting expressions for the ESP at large 



J > fi , 



w j and 



small J < fi 



,u>y [Wj: is the correla- 



tion rate of V(t) fluctuations (see below)] turn out to be 
similar to those obtained above in the cage model>iZr— 



V = V s = (w Pf ) 



for J > fio,Wj, (3.17) 



V = V.,„ 



1 + ujf 

J \w f ){l +ojiy 1 



Here 



fi /tty and Pf 



{\y±\ 2 )\(^_ 

\aD p 



(3.19) 



with Wf — w^, ( (i = H, Z), being the decay rate of the 
correlation function (V^tyV^O)) (defined in Sec. IIIA). 
Similarly to the cage model p R — l s /r i « l s /d is 

the probability of (diffusion controlled) reaction in S- 
state with the rate k r , given by eq. (|2.8I) . in which l s = 
d(K r d/ D p )[l + (n r d/ Dp)}^ 1 is the reaction radius. 



in which ip u = arctan(w_j_) + (2cu c ) 1 ln(l + uj 2 .) with 



J± =(fi ±2J )/u; c , 
uj c = n /w c and p c 



, ((\V ± \ 2 )\(^d 2 



w c w d 



(3.14) 



Here p Rc — n R /(w d + k r ) is the probability of reaction in 
5-state and V c is the volume of the cage. 

Note that in the presence of V(£)-fluctuations the rate 
formulas (|3.13p and (|3.14p are still valid but w c — w d + 
^k r + io„, ( p = H, Z), where is the rate of decay of 
correlation function (^(4)^(0)), introduced above. 



C. Some general properties of D — D ESP 

The important specific features of the ESP for differ- 
ent limits and mechanisms (discussed above) are conve- 
niently represented in terms of the dimensionless func- 
tions $*(o;), in which the superscript p specifies the 
ESP mechanism (see Sec. IIIA): HFI (jj, = H) and AZI 
(p = Z), while the subscript \ indicates the value of Jo: 
large (x = s) and small (x = w). In accordance with 
above-obtained expressions, one can write 



(3.20) 
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Here and w 7 are denned in eqs. (|3.14[) and (|3.19p 

for the cage (7 = c) and free diffusion (7 = /) models, 
respectively. The parameters (v 2 ) are the dimensionless 
squares of interactions defined by formulas 

= ^? {B2 i^f b \ {<) - g ° : ^/» :g » , (3-2i) 

with w , (7 = c, /), given in eqs. (|3.14l) and (|3.19[) . 

a. Strong exchange interaction limit (x — s). In the 
limit of strong exchange interactions (large J ) we get for 
the HFI (p = H) and AZI (/.i = Z) mechanisms 

$£(w)=w/(l + w 2 ), $;H= W 3 /(1+^ 2 ), (3.22) 

where the parameter P s is written for the cage (7 = c) 
and free diffusion (7 = /) models as follows: 

P s = {^){w c /w d ){^d 2 /aV c ) for 7 = c, (3.23) 
P s = (±n)(dw f /aD p ) for 7 = /. (3.24) 

b. Weak exchange interaction limit (x — w). In the 
weak exchange interaction limit (small J ) one obtains 

$-H= W /(l+ W 2 ) 2 , <I>-H=c 3 /(1+^ 2 ) 2 , (3-25) 
with 

P w = (8/tt)(J /w^)P s , (3.26) 

both in the cage and free diffusion models (for 7 = 0,/). 

The functions $^(w) are shown in Fig. 2. Notice that 
all dependences except $J(w), are non-monotonic 

with maximum at oj ~ 1. 

Note that in the free diffusion model the formulas for 
the ESP are derived under the assumption of not very 
high mobilities of particles (fairly small diffusion coeffi- 
cients) so that £ n = flod 2 /D > 1 J£ In the opposite limit 
£ Q = flod 2 /D p < 1 the ESP is essentially determined by 
the parameter £ w = w c d 2 /D p . For £ w > 1 we get the 
results presented above independently of the value of £n . 
For £ w -C 1, however, the behavior of at £ w < I 

becomes non-analytic: for example, ~ w _3 / 2 J^ 

IV. T-D QUENCHING 

The T—D quenching is another example of spin selec- 
tive processes, leading to the generation of the ESP of 
D-particles. The kinetic scheme of the geminate T — D 
process can be represented in the form (|2.I[) with X a = T 
and Xjj = D (in OSCs D-particle is associated with 
trapped electron). 

The quenching is believed 2 ^ to result in the transition 
to X R — (S D*)d, where S is the molecule in the ground 
(singlet) state, and D* is the excited doublet (free elec- 
tron), i.e. the quenching rate is expected to be non-zero 
only in the doublet (|-D ±1/2 )) states of T — D pair, cor- 
responding to the total spin S — 1/2. In this case the 




FIG. 2. The functions <&*(w), which determine the ESP de- 
pendence on the Zeeman splitting cj = fio / w ~i in cage (7 = c) 
and free diffusion (7 — c) models [see eqs. (13 . 22 (I and (|3,25[) ]: 
(1) (2) %(u), (3) and (4) [see eqs. 

(|3.22[) and (|3.25f) ] . In these functions the superscript \ and 
subscript p specify the exchange interaction strength and the 
mechanism of the ESP generation, respectively, (x = s,w and 
a = H,Z). 



quenching operator K R is represented in the form (|2.3p 
with P R =P D = \D +1/2 )(D +1/2 \ + \D_ 1/2 )(D_ 1/2 \. 

It is clear from the above comments to the general 
scheme (|2.1I) that, in principle, there are two types of 
-D-particles involved in the process under study: D and 
D* [escaped from the reaction (D) and resulting from 
reaction (D*)]. In OSCs these particles are associated 
with trapped electrons and electrons in in the conducting 
band, respectively. 28 The T—D quenching gives rise to the 
ESP of both particles. In the analysis, however, we will 
assume that the kinetic and MFE parameters of these 
types of -D-particles are different, i.e. the contributions 
of them to the ESP are experimentally distinguishable. 

In our work we will concentrate on the evaluation of 
the ESP of D* particles. 

The methods of estimating ESP of both D particles 
for liquid phase processes, i.e. in the diffusion model of 
relative motion of particles, have already been discussed 
in a number of papers i 17 ' 20 ' 21 Below we will analyze some 
specific features of the ESP generation as applied to dis- 
ordered semiconductors, considering more thoroughly the 
cage model, as an example. 
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A. Interactions 

In our analysis we assume that the spin Hamiltonian of 
the T—D pair (similar to that for D—D systems) can be 
represented as a sum H(r, t) — H (r)+V(t), where H (r) 
and V(t) are the steady state and randomly distributed 
(and may be fluctuating) parts, respectively. 

The steady state part 



H (r) = H z + H ex (r), 



(4.1) 



diagonal in the basis of the total electron spin S = S T -f 
S D , is a sum of the Zeeman and exchange interactions: 



H z = lj S z and H ex = J(r)(| + 2S T S £ 



(4.2) 



in which S z — S T + S D , loq — gfiB (we assume that 
9d = 9t =5 = 2), and J(r) is given by eq. (|3.5|) . The 
terms U M (r), (M = Q m ,D m ), of the Hamiltonian H (r) 
are displayed in Fig. 3. 

Transitions in the regions of crossing of U M (r)-terms 
are determined by the fluctuating zero-field-splitting 
(ZFS) interaction 



V = V T =D t {S 2 t , -lS 2 )+E T (S 2 



s 2 



(4.3) 



where S T . , (j = x',y', z'), are the projections of the spin 
of T-exciton on the eigenaxes of the ZFS tensor. Usually 
D T 3> E T , that is why below we will neglect the terms 
~ E T , taking V T ~ D T S 2 t . The fluctuations result from 
hopping motion of T exciton over sites with randomly 
oriented directions of axes x',y',z'. These fluctuations 
show themselves in those of interaction matrix elements 
V T , (f) = (M\V T (t)\M'), correlation function of which 
is assumed to be of the formSZ Or* (i)K (0)) o , = 

V T MM' K ' MM' " n ' 

(\V 2 I ) e _w r* (here M,M' = Q ,D and 

MM' 

averaging over orientations of the ZFS eigenaxes) 



) , denotes 



B. ESP generation kinetics 

In T — D system the ESP is assumed to result from 
quantum transitions in the regions of crossing of terms 
Un (r) and Un (r) with reactive terms Ud (r). 

V ±3/2 V ' V ±l/2 V ' ±1/2 V ' 

These regions are indicated in Fig. 3. 

Similarly to the D—D ESP, the ESP in T—D quenching 
will be evaluated for the equilibrium initial spin state 



Po = Pe = \{P D + p Q ), 



(4.4) 



where P D and P are the operators of projection onto 
the subspaces of D- and Q-states, respectively: P = 

£ m=±I/ J£LXAJ and p q = E_ ±1/2 , ±3/2 IQJ<QJ- 

The ESP for other initial spin state are, actually, pro- 
portional to that for p E i2£ 

In the majority of T—D processes the above-mentioned 
ZFS-induced quantum transitions can be calculated per- 
turbatively in the second order in V . In this approxima- 
tion the effect of transitions is represented as a sum of 




FIG. 3. The schematic picture of terms U M (r), (M = 
D ±i/2 1 Q±i/ 2 ' Q±3/2 )' °f the s P in Hamiltonian H [r) [eq. (|4.1jl ] 
of the T — D pair (for J > 0), in which r tl 2 and r t3 are the 
coordinates of the regions of efficient quantum transitions. 



contributions of all pairs of coupled terms, which can be 
evaluated with formulas derived above for the D—D ESP 
PJJL It is worth noting that these formulas represent the 
effect of quantum transitions on the D—D ESP of escaped 
articles. However, similar expressions are, certainly, valid 
for the product particles as well^ 

Analysis of quantum transitions between terms (dis- 
played in Fig. 3) shows that the contributions to the 
ESP V* of the D*-particle can be represented as a sum 
of three contributions of type of (13.121) : 



V = P 



2 sv 

3 I 1/2,3/2 



SV 



1/2,1/2 



5V 



1/2,-3/2 | ' 



in which SV , are defined by 



SV 



V 



V 



[with to = 1/2, ml = 1/2, 3/2], where V D 
contribution of the transition D <— Q , . 



(4.5) 



(4.6) 
is the 



1. Evaluation of ESP in the cage model 

In the cage model we can obtain the expression for 
SV , in terms of the distance dependent solution 
p I (r) similar to that derived above [see eq. (|3. 121) ] 



SV mr 

In this formula 



w 2 + Q?(r) w 2 c + ttl(r) 



(4.7) 



(4.8) 



where the parameter p = ^(K R w c )/[w d (w d + k r )] is iden- 
tical to that in eq. (|3.12[) . and 



0^7-)= a n Q ±2J(r) 



(4.9) 
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With Q ? ±1/2 ,±3/2 Q! ±l/2,^l/2 

r2 



1 and a ±1/2 T3/2 = 2. 



/2,±3/2 I' 

4 n2 



As for (\V 2 J), for m and m', corresponding to the 
contributions (|4.6|) . they are written as (|V^ 2 

1E D 1 = lW> ^d (|y ± 2 1/2 , T3/2 |) = ^ 

Note that similar to the D—D process in the presence 
of ^-fluctuations the dephasing rate is written as w c = 
Wd + 5 «k + W T [rather than as that defined in eq. p.20jl ] , 
where w T is the decay rate of the correlation function 
(V* (f)V T (0)) o; , defined in Sec IVA. 

T MM> J AfAf' V " n 

Substitution of all terms (|4.8[) into eq. (|4.5p yields the 
final expression for the ESP of D* (see below). 



2. Combined formula for ESP, predicted by two models 



The T — D ESP, predicted in the above- analyzed cage 
model and previously studied free diffusion model ) 17 ' 20 i 21 
can be represented in one combined expression similar to 
that (f3T20]> obtained for the D-D ESP 

n = lv Ri (vl)P x [^ 1 ) + 2^{2u Jl )], ( X = s,w). (4.10) 

In this expression (u 2 ) = j|(Z3 2 /u> 2 ), the functions 
$*(w 7 ) are defined in eqs. (f3^2|) and (pT25|) . and the 
parameters P are given in eqs. (|3.23|) . (|3.24p . and (|3.26|) . 

The difference between the cage and free diffusion ESP 
lies in the parameters p^ and oj^: for the cage ESP (7 = 



c) P Rc = K R /(w d + k r ) and uj c = Vl Q /w c [see eq. (|3.14|)], 
while for free diffusion ESP (7 = /) p = (k d/D p )/[l + 
(K r d/D p )} and ajy = f2 /uy with ity = w T [eq. (|3.19|) ] . 

It is worth noting that the dependence of V t on uj y , 
(i.e. on f2 ), predicted by both the cage and free diffu- 
sion models, is determined by that of i.e. function 
V^Wy) is non-monotonic the maximum of 7^(u!^) is lo- 



cated at ui. 



1 (see Fig. 2) 



DISCUSSION 



The obtained formulas allow us to analyze easily the 
specific features of the ESP in D—D and T—D processes. 



A. General remarks 

Before estimating the ESP some remarks are needed 
on the validity of models and approximations applied. 

1. So far in our work we have discussed the absolute 
value V of the ESP. As for the ESP sign, it can be charac- 
terized by the sign factor S v = ±1, defined by the relation 
2(S Dz ) = S V V, in which the average (...) is made over 
spin state and distances (see Sec. II). The factor S v can 
easily be obtained from the qualitative analysis of quan- 
tum transitions in D — D and T—D systems. In general, 
Sp is essentially determined by the sign of the exchange 
interaction S v ~ sign(J ). Below we will briefly discuss 



the ESP sign in the most realistic case Jo > (for Jo > 
the terms are shown in Figs. 1 and 3). 

a) In the D—D process the ESP of survived D-particles 
is determined by STl-transitions, leading to the decrease 
of the population of T_-state, i.e. to S v = +1. 

b) In the T - D quenching the ESP of D* particles 
results from quantum transitions, associated with contri- 
butions (|4.6|) resulting in larger population of D +1 , 2 -state 
(because of stronger coupling between initial Q states and 
this D-state). Therefore for D* particles S T = +1 as well. 

2. In the proposed theory the ESP V is evaluated in the 
second order in £ 7 = (\V 2 \)/w 1 <C 1, (7 = c, /), which 
implies that V <C 1. This does not mean, however, that 
the generated ESP is small, in general. For example, the 
effect of attractive interparticle interactions can lead to 
strong increase of the ESP (see Sec. V.B.2). 

3. In the cage model the spin selective reaction rate 
k r is suggested to be independent of the distance r. This 
assumption is, in principle, not quite consistent taking 
into account that J(r) is treated as distance dependent. 
Notice, that the dependence k r on r could be analyzed 
within generalized formulas [see eqs. (|3.12[) and (|3.12[) ]. 
though, the modification of formulas would result in the 
strong complication of the expression for V . To justify 
the simplifying assumption on k r independence of r it is 
worth noting that the manifestation of K H (r)-dependence 
in V is fairly weak. It is especially weak in two above- 
considered limits of strong and weak exchange interac- 
tion, most important for further applications. 

a) For strong exchange interactions the quantum tran- 
sitions, determining the ESP, are highly localized in the 
regions of term crossing (Figs 1, 3). The distance depen- 
dence k r (r) in these regions can be neglected. 

b) In the weak exchange interaction limit the depen- 
dence of V on the rate k r results from that of the proba- 
bility of reaction in the cage p R ~ k r , (7 = c, /). In the 
free diffusion model the weak effect of the r-dependence 
of reactivity is demonstrated in earlier worksJi~— Here 
we will discuss the cage model, in which this effect shows 
itself in the average over r in formula (|3.11[) for V, i.e. 
in the average {p R J) r , where p R = K R /(w d + k r ). It is 
seen that in (p Rc J) r the K R (r) dependence manifests itself 
only weakly. For example, in the low reactivity limit, 
when p R w K R /w d and the manifestation is strongest, 
the dependence K R (r) ~ e~ ar , identical to that of J(r) 
[see eq. (|3.5[) ]. results only in two-fold reduction of the 
ESP. This effect is too small to be of any importance for 
our semiquantitative study of the problem. 



B. Estimations of the ESP 

1. Bulk processes 

In our analysis we have considered the ESP generation 
in geminate processes. It can be known, however, that 
the geminate ESP V, calculated for the equilibrium ini- 
tial spin state, is directly related to the rate JC of ESP 
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generation in the corresponding bulk processes^ 



/C = k V with fc = 4irr c D p 



(5.1) 



(note that for T—D system V = V^). In this formula D p 
is the coefficient of relative diffusion and r c is the radius 
of the cage (r c ~ A) in the cage model, or the radius of 
efficient quantum transitions (r c ~ d) in the free diffusion 
model. Our assumption that d ~ A (Sec. II) means that 
in both models the estimations for r c are close. 

Note that the parameter r c , whose meaning in the dif- 
fusion theory is clear; 17 ' 18 ' 20 ' 21 is not quite well defined 
in the cage model. This uncertainty of the cage-model 
definition of r c is, however, not very strong and is not 
essential for our semiquantitative analysis. 

Sometimes it is convenient to relate the rate of ESP 
generation to the reaction rate: 



(PsP^K, (7 = c,/), 



(5.2) 



where p^ is the probability of reaction in the reactive 
state [see eqs. (|3.14[) and (|3.19[) ] and p s is statistical 
factor, representing the relative number of reactive spin 
states: for D — D and T — D pairs we get p s = 1/4 and 
p 3 = 1/3, respectively. In eq. (|5.2p we neglected small 
contribution of quantum transitions to the reaction rate. 

Formulas (|5 . 1 1) and (I5.2f) are fairly suitable for studying 
the efficiency of the ESP generation in different types of 
bulk processes. The explicit expressions for the ESP are, 
naturally, determined by the particular kinetic scheme of 
the process, and should be discussed as applied to the 
interpretation of specific experiments. 



2. Effect of attractive interparticle interaction 

The proposed theory is based on the assumption that 
particles X a and X a [see eq. (|2.ip ] undergo free stochas- 
tic motion. In reality, however, the ESP can be strongly 
affected by the interparticle interaction U(r) = k B Tu(r). 

Especially strong effect of the interaction is expected 
in electron-hole (e — h) recombination, kinetics of which 
is known to be markedly influenced by the Coulomb e — h 
interaction. Significant interaction effect can also be ob- 
served in some T—D pairs, in which D-particle is the elec- 
tron (e) or hole For these pairs U{r) ~ — ^Oj, i? 2 (r), 
where Oj, is a polarizability of T state and E(r) — e/ (er 2 ) 
is the strength of the Coulomb field of e (or h) particle 
in the semiconductor with dielectric constant e. 

In the presence of the well-type attractive potential 
U(r) the ESP V c can approximately be estimated within 
the simple kinetic scheme, taking into account the ef- 
fect of the cage within the well , 24 ' 25 which hereafter is 
denoted as [/-cage. The size R c of the [7-cage, deter- 
mined by the relation \U(R C )\ = k B T, is expected to be 
significantly larger than the size r c of the primary cage, 
considered above in the cage model: R C ^S> r c ~ A. 

The kinetic scheme of the process in the [/-cage (in 
the well) can also be represented as (|2.1j) . but with the 



reaction (W r ) and dissociation (W d ) rate a 24 ' 25 
W r = K r /V and W d = {KjV u )e- u », 



(5.3) 



in which u a = U a /{k B T) is activation energy of dissoci- 
ation (the well depth), K = AnR c D p is the bimolecular 
rate of capture into the well, and V u is the volume of the 
[/-cage. The volume can be estimated as V v « (47r/3)i?^, 
where iJL is the radius of the thermal distribution within 
the well, i.e. U{R U ) « U b + k g T, with U b = -U a being 
the energy of the bottom of the well. 

In the [/-cage the ESP generation is described by the 
rate 



W c = TW , where W = k /V v 



(5.4) 



The kinetic scheme (|2.1|) predicts exponential kinetics 
of the decay of [/-cage population, but with different de- 
cay rates for different spin states of pairs. In particular, 
for non-reactive states (T and Q states for D — D and 
T — D pairs, respectively), whose depopulation kinetics 
determine the ESP generation, the decay rate is Wd- In 
such a case the [/-cage affected ESP V c is written as 



r c = r(w /w d ) = v(r c /R c y 



(5.5) 



(recall that for T — D system V = V^j. Similarly the U- 
cage affected ESP generation rate K. c is given by JC C — 
V C K = K,e Ua . As for the decay rate of reactive states, it 
is equal to W r + Wd and, therefore, the [/-cage affected 
reaction rate is represented as K. Tc = [W r / (W r + W d )]K c . 

It is seen from eq. (|5.5p that the attractive interparticle 
interaction can result in the strong increase of the ESP. 

Noteworthy is that formulas (|5.5[) are obtained for not 
very deep wells, when r\ — W c /W d < 1 and one can ne- 
glect the contribution of quantum transitions to the re- 
action rate W r . To find V c and K. c in the opposite limit 
t] > 1, one has to solve complicated kinetic equations, in 
general. Some simple estimations can, nevertheless, be 
made in the realistic case of strong J(r), large Zeeman 
splitting w > 1, and high reactivity (when p R w 1). In 
this case only quantum transitions between non-reactive 
and reactive states, localized in crossing regions, con- 
tribute the ESP. These contributions can be easily eval- 
uated for both D — D and T — D reactions: 

a) For D—D system only T_ —> S transitions contribute 
to the ESP so that V r can be written as 



|w ST /(W d + W ST ), 



(5.6) 



where W ST _ = W S ^ T _ = 4PW is the transition rate, 
whose value is obtained by comparing the expression 
([53)1 with eq. (153)) in the limit W ST _ < W d . 
b) For T — D system similar analysis leads to formula 



V = V 



l[W 2 /(W d + W 2 ) + W_/(W d + W + )} 



= (nw )/[Wd + (9/2m>v ] 



(5.7) 



in which W ± = W 3 ± W t , with W 1 



|^W and W 3 



W r 



= W D ^ Q 
2W 1; and W 2 
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W n „ = 3Wi. Similarly to the case of D-D 

system the expression of the rates Wj, (j = 1 — 3), in 
terms of the ESP V = V* is found by comparison of eq. 
(I5J1) with eq. dHU) (for V c = 7> J at }% < W d . 

Formulas (15.61) and (15. 7[) . looking quite natural in the 
diffusion model, need some comments as applied to the 
cage model. These expression imply that within the pri- 
mary cage (at r < r c ) a large variety of relative positions 
r of particles, including those in the term-crossing re- 
gion, are accessible. Note that possible migration of par- 
ticles over these positions can result in ^(^-fluctuations 
and, thus, to some change of w c (see Sec. IIIB.l and 
IVB.l). However, in the considered case of strong J(r) 
and large w 7 this change does not lead to any change 
of the ESP32&2 1 . Noteworthy is also that the above- 
mentioned large variety of r in primary cages can be re- 
alized as a result of reencounters of reacting particles at 
different points in the [/-cage. In such a case the ESP is 
quite properly described by the original eqs. 1)5. 6p . (|5.T[) . 

3. Magnitude of the ESP 

As we have already noted above, the dependence of the 
D-D ESP V{uj^) ~ 3>X(w 7 ) on the dimensionless Zee- 
man splitting oj 1 = Cl /w^, (7 = c, /), is non-monotonic 
for both HFI and AZI mechanisms of ESP generation 
(see Sec. IIIA) and both models of relative motion: the 
cage (7 = c) and free diffusion (7 = /) [except the case 
fi = Z, x = s (|3.22l) ]. Similar non-monotonic behavior 
of V^ljS) is observed in T — Z?-quenching [eq. (|4.10|) ]. 
In all these cases the maximum value V m = V(ui„ m ) [or 
V m = "P*(w 7m ) for T — D process] is located at u>„ m ~ 1. 

Such a behavior of V(ui^) allows one to simplify the 
analysis of the ESP by restricting the discussion to max- 
imal values V m of the effect. In what follows, for cer- 
tainty, we will estimate V at o> = 1, for which 3>£(1) ~ 
[$*(l) + 2$*(2)] ~ 1, so that the maximal values V m for 
D — D and T — D processes can be written as 

~ PaJ$)Px, (M = H,Z,T and X = s, «>), (5.8) 

In this formula p^ is the probability of reaction in reac- 
tive spin states [see eqs. (|3.14p and (|3.19[) ]. The param- 
eter (v£) is the dimensionless coupling of terms of H (r) 
[eqs. (|3~2T|) and (j¥7TU|) . As for P x , it describes the de- 
pendence of the ESP on the strength of J(r) [\ — s an d 
X — w for strong and weak J, respectively). Formulas 
for P x , different in two considered models of motion, are 
presented in eqs. (1535)1 , 1)5341) and 1)535)1 . 

Of certain interest is the above-mentioned case of AZI 
mechanism (fx = Z) and strong J(r) (x = s), in which 
$| (oj) is a monotonically increasing function. Below it 
will be analyzed more thoroughly. Here we only note that 
in this case at large lu ~ B the ESP can be very large. 
Now we will estimate the parameters in eq. (|5.8|) : 
1. The largest ESP is, naturally, expected in the case of 
strong reactivity resulting in a high reaction probability 



p^ . In what follows we will assume that p n ~ 1 . 

2. The parameter (v 2 ) (fx = H,Z,T) is suggested to 
be small: (v£) < 1, to ensure the validity of the applied 
perturbation approximation (in V). Noteworthy is that 
for AZI mechanism (fx = Z) the validity of this approxi- 
mation implies also not very large value of the magnetic 
field B (to be sure that the value of (V^)/w 2 is small). 

3. The parameter P x [as it follows from eqs. (|3.23l) . 
(I3.24[) and ()3.26p ] is smaller than 1, but is not very small. 
The fact is that, in reality, w e ~ w d ~ uy and D p ~ 
X 2 /w d ~ d 2 /w d , so that the small value of P x results 
from the relation P x ~ (ad) -1 , in which ad < 10. 

Combining the estimations of these three parameters 
we can obtain the ESP, but to calculate it more accu- 
rately we need to specify the process under study. 

For example, let us consider the cage model prediction 
for the ESP V = V* of D* particles in T — D quenching. 
The value of the parameter (v£) for the T — D system is 
determined by the ZFS constant D T and w c ~ w h . For 
realistic values D T = 3 ■ 10 9 s _1 and w h = 10 9 s _1 one 
obtains (v 2 ) = j^(D 2 /w 2 ) ~ 0.4. Taking into account 
the above estimations p^ ~ 1 and (ad)" 1 > 0.1 we get 
(for u> c = 1) the value 7\ > 4 • 10 -2 , much larger than 
the thermal ESP V„ th ~ 1.4- 10~ 3 (at room temperature). 

This value of is obtained in the case of strong 
exchange interaction ( x — s), predicting largest ESP. Of 
course, in the opposite limit (x = w) the ESP is smaller 
with the factor J /w -C 1 [see eq. (|3.26l) ] . 

The attractive T—D interaction can result in the strong 
increase of the [/-cage affected ESP 7^ : V <fa ~ V 3f e Ua , 
where u a is the depth of the potential well ([/-cage), as 
it follows from eq. (|5.5p . In particular, even in the case of 
fairly shallow attractive well of the interaction potential, 
corresponding to u a = 2.0 (fcgT), this interaction gives 
rise to the significant (about one order of magnitude) 
increase of the ESP V* ■ Taking into account the above- 
obtained estimation of we expect, that the interaction 
affected T—D ESP can be fairly large: V tte ~ 0.1. 

Similar estimations can also be made for the D — D 
ESP. Because of this similarity we are not going to do 
them here, but restrict ourselves to some comments on 
the AZI mechanism of the ESP (ll = Z). The fact is that 
values of (v 2 ) [see eq. (|3.2ip ] typical for organic semicon- 
ductors are very small: (v 2 ) < 10~ 5 [due to small g':g' 
1)3.8)1 ]. and for not very strong magnetic fields B, corre- 
sponding to uj 7 < 10 3 , the ESP is fairly small: V < 10~ 3 . 
It can be significantly larger, however, in the presence 
of the attractive interaction ([/-cage), as it follows from 
above relations. Note also, that in some organic semi- 
conductors, doped with heavy atoms, the Zeeman inter- 
action can be strongly anisotropic^ resulting in fairly 
large (v 2 ) > 10 -4 and, thus, in the ESP much larger 
than the thermal one. 

Concluding our analysis of the [/-cage effect on the 
ESP, we recall that the applied formulas 1)5. 3p and 
1)5. 5|) are valid for relatively shallow wells, when 77 = 
(■PW )/>Vd < 1. For deeper wells, for which 77 > 1, a 
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fairly large £/-cage effect is expected, though the calcu- 
lation of the ESP in this limit is a rather complicated 
problem, in general. Fortunately, in the realistic case of 
strong J(r) and large Zeeman splitting to it is simplified, 
as has been shown above [eqs. (|5.6|) and (15. 7[) ]. Obtained 
formulas predict that in the limit r\ 3> 1 the ESP is very 
large: V c (rj > 1) = 1/4 and P^(?] > 1) = 2/9 for D-D 
and T—D processes, respectively. 

Naturally, all above discussed effects manifest them- 
selves in the rates of ESP generation in corresponding 
bulk processes as well. We are not going to analyze this 
manifestation. Such an analysis would require the specifi- 
cation of the process under study, which should, perhaps, 
be done as applied to particular experiments. 

Here we will only restrict ourselves to a general re- 
mark on the relation between geminate and bulk ESP. 
Note, that the ESP can be analyzed in terms of the ra- 
tio of the ESP generation QC) and reaction (JC r ) rates: 
Rf- = fC/K, r = (PsPj,) r. Similar relation is valid in the 
presence of the attractive potential [see eqs. (|5.5[) ]. All 
these relations show that the analysis of the ESP genera- 
tion in bulk processes provides the information identical 
to that obtained by studying the geminate ESP. 

VI. CONCLUSIONS 

In this work we have analyzed the specific features of 
the ESP of charge carriers (electrons and holes) in D—D 
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